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Abstract. The mean- field pictures based on the standard time-dependent variational 
approach have been widely used in the study of nonlinear many-boson systems such 
as the Bose-Hubbard model. The mean-held schemes relevant to Gutzwillcr-like trial 
states \F), number-preserving states |£) and Glauber- like trial states \Z) are compared 
to evidence the specific properties of such schemes. After deriving the Hamiltonian 
picture relevant to \Z) from that based on \F), the latter is shown to exhibit a Poisson 
algebra equipped with a Weyl-Heisenberg subalgebra which preludes to the |Z)-based 
picture. Then states \Z) are shown to be a superposition of A/"-boson states |£) and the 
similarities/differences of the |Z)-based and |£)-based pictures are discussed. Finally, 
after proving that the simple, symmetric state |£) indeed corresponds to a SU(M) 
coherent state, a dual version of states \Z) and |£) in terms of momentum-mode 
operators is discussed together with some applications. 



PACS numbers: 03.75.Fd, 03.65. Sq, 03.75Kk 
1. Introduction 

The semiclassical formulation of many-mode boson models based on the coherent-state 
(CS) method pQ, [2] has proven to be an effective tool to describe the behavior of 
interacting bosons in many situations [3J-JH]- Such models, usually represented by 
a second-quantized Hamiltonian in terms of boson operators with 
standard commutators [a m , af] = 5 m i, exhibit a dramatic complexity owing to their 
many-body nonlinear character. The combination of the CS method with the application 
of standard variational schemes allows one to circumvent this problem by reformulating 
model Hamiltonians into a mean- field (M-J 7 ) picture [B], [12] in which the Schrodinger 
problem for variational trial states |<&) = \<f)i,(p2 ■■■) describing the system quantum 
state is reduced to a set of Hamilton equations governing the evolution of parameters 
A very standard choice [13J is |$) = \Z) = Yli \ z i) where, for each mode aj, 
state \zi) is a Glauber CS satisfying the defining equation a m \z m ) = z m \z m ) and (f)j 
identify with CS parameters Zj = ($|aj|$). Similar schemes have been developed for 
magnetic and multi-level atomic systems (6j, [H] where |$) is a product of spin CS. The 
important feature is that dynamical variables <pj are, at the same time, the expectation 
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values of the Hamiltonian operators and thus are providing information on microscopic 
physical processes. Significant examples are found within the physics of ultracold bosons 
confined in optical lattices where this M.T formulation successfully describes complex 
dynamical behaviors [15]-[25]. There, CS parameters usually represent either the on- 
site superfluid order parameters of local condensates or the expectation values of local 
operators such as a\a m and giving informations on space correlations and boson 
populations, respectively. 

The most general version of this M.T picture, however, is achieved by using a 
Gutzwiller-like [26J trial state |<&) exhibiting yet a factorized semiclassical form, but 
involving constituent states more structured than CS. These are 

M M oo 

i^=ni F <>=nE&i*>*> c 1 ) 

i=l i=l rt;=0 

where \Fj) has replaced local state \zi) in \Z), rii\n)i = n\n)i and M is the boson- 
mode number namely, for many models, the lattice-site number. This choice ensures an 
improved description of microscopic processes in the sense that, for each mode, infinitely- 
many variational parameters are available now in place of the M parameters Z{ of \Z). 
The |F)-based approach has been applied in studying the dynamics of Bose- Hubbard 
(BH) model [27], [28] as well as its zero-temperature critical properties [29j-[32j. 

Recently, a third variational scheme has been considered [23] , [25J to approach the 
dynamics of many-mode boson models, where state |$) is assumed to have the form 

|Ar,O = (AT!)- 1 / 2 (A + )^|0,0,...,0), (2) 

where A + = Y2iLi£,i a t an d the constraint Y2fLi |£«| 2 = 1 ensuring its normalization. 
Different from states \Z) and \F), the distinctive property [33] of states |£) (A/* will be 
often implied in |JV, £)) is to diagonalize, by construction, the boson-number operator 
N = £\ Hi whose eigenvalue can be easily shown to identify with index M. Hence, states 
|£) having M as a good quantum number naturally embody the property [N, H] = 
characterizing usually many-mode boson Hamiltonians H . This valuable feature reflects 
in turn an even important fact that states [2] actually coincide with the CS of group 
SU(M) where eigenvalue M is the index labeling the representation of SU(M). 

The structure of formula [2} however, appears quite different from the (standard) 
group-theoretic form of SU(M) coherent states. The standard definition [I], in fact, 
states that |£) = g\Vt) with g G SU(M) where \fl) is an appropriate extremal state. 
Then a CS is generated through the exponential action of an algebra element a G su(M) 
such that g = exp(ia) where a is in general a linear combination of su(M) generators. 
A well-known M boson-mode realization [2] of SU(M) CS is, for example, 

ff |n> = r(c)|n>, |o> = |o, ...,tv,o, ...> 

where T(£) = e ta is the displacement operator, a = X^ m (C a m a ^ + a / a m), Ct £ C, 
and ne\Q) = 5£ m Af\Q). Such a definition has been (and is) currently in use to 
represent quantum dynamical processes in microscopic systems within Quantum Optics, 
Condensed-Matter theory, and Nuclear Physics (see [2], [31] and references therein). 
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Except for the case M = 2, where states [2] are easily related [2 J to definition g\Q) 
{ see | Appendix A[), for M > 3 the connection of formula [2] with the group-theoretic form 



of SU(M) CS is less direct owing to the difficulty in disentangling the group element g 
through the Baker- Campbell-Hausdorff decomposition [2]. As noted in [35], where this 
issue has been investigated, even if SU(2) CS have been widely used in the literature, 
not much work has been devoted to find realizations of SU(M) CS more practicable in 
physical applications. In this respect, indeed definition [2] has proven to supply a useful 
tool in applications. However, despite states [2] being tacitly presented as CS of group 
SU(M) in various papers, their connection with the group-theoretic definition within 
CS Theory is far from being evident. We thus devote some attention to this particular 
aspect even if it might be known to authors [35]- [38], [39] involved on mathematical 
aspects of CS. 

This paper is aimed at comparing the M. T schemes based on states \F), \Z) and 
widely used in applications to boson systems. We enlight some formal aspects concerning 
both the implementation of the (time-dependent) variational approach within such 
schemes, and the representation of trial states in terms of CS. We emphasize that 
some parts of our discussion have a review character and involve well-known theoretic 
tools. Nevertheless, a direct comparison among these three M. T schemes has never been 
presented to our knowledge. We feel that such a comparison can elucidate their specific 
properties and advantages at the operational level prompting as well their applications. 

A first objective of this paper is to evidence how the variational schemes based on 
state \F) and \Z), respectively, are related to each other. In section El after reviewing 
the variational procedure that amounts to reformulating quantum Hamiltonian models 
in terms of effective Hamiltonians and of the corresponding dynamical equations, we 
focus on the formal derivation of the Hamiltonian picture relevant to \Z) from the one 
based on \F). We show that the Poisson algebra of variables f^, is naturally equipped 
with a (classical) Weyl-Heisenberg sub-algebra preluding to the |Z)-based picture. 

A second objective is to relate SU(M) CS |£) to states \Z) and \F). In section [3j 
following [13] . we show that Glauber- like trial state 

M co n 

\Z) = TT \zi) , \zi) = e^-^|0> 4 = e-l^l 2 / 2 V ^j=\n)i , z t e C 

can be expressed as a superposition of SU(M) states |£) thus making evident why the 
dynamical equations obtained in the |Z)-based scheme have essentially the same form of 
those obtained in the |£)-based scheme. To illustrate this situation we derive the M.T 
dynamics relevant to the BH model showing (see Appendix D ) how the use of form [2] in 
place of the standard SU(M) CS definition is extremely advantageous. In this section 
we also display an explicit way to relate state [2] with the standard form |£) = g\Q) 
of CS theory involving one among M possible (equivalent) choices of extremal vector 
| Q) and of the relevant maximal isotropy algebra. Finally, in section [H we discuss 
the property of "duality" inherent in space-like states \Z) and |£) defined in terms of 
ambient-lattice operators af showing how they can be easily rewritten as momentum-like 



Some remarks on the coherent- state variational 



4 



states involving momentum modes bt . We exploit this property to construct Schrodinger 
cat states in terms of states |£) and to show that they exhibit specific momentum 
features. 



2. Mean-field approaches based on states \F) and \Z) 

In order to compare the |Z)-based approach with the |F)-based approach we refer to 
the M.T dynamical equations stemmed from such schemes for the well-known Bose- 
Hubbard Hamiltonian [6]-[Tn]. Optical-lattice confinement shows that real boson lattice 
systems are effectively described within the Bose-Hubbard picture. The corresponding 
Hamiltonian HQ], defined on an M-site lattice, is 



H = ^ J^(nJ - rij) - T H a t a t ( 3 ) 

3 <j£> 

where = a\ai (i — 1, M) and a h a\ obey the standard commutators [a m , af] = 5 mi . 
In the hopping term Yli<ji> = \ J2j S^ej' wnere ^ labels the nearest neighbour sites of j, 
and Tj£ = T#j. This model well represents the boson tunneling among the potential wells 
forming d-dimensional arrays (d = 1, 2, 3) through the hopping amplitude Tje, and takes 
into account boson-boson interactions by means of parameter U. For one-dimensional 
homogeneous arrays the hopping term reduces to T£j(a+ +1 a,- + a J+ ia+). 

The M.T dynamics relevant to trial state [1] is easily derived. State \F) = Ui\Fi), 
where = fn\n)i and |n)» is such that a^\n)i = \/l + n\n+l)i, a^nji = y/n\n—l)i, 
obeys the normalization condition (F\F) = IT [Y^ n \fn\ 2 } = 1- The application of the 
time-dependent variational principle [2], [6] amounts to deriving dynamical equation for 
parameters f l n by stationarizing the weak form of Schrodinger equation (\I/|S|\I/) = 
where S := ihdt — H. In order to illustrate this procedure, we write explicitly the 
expectation value of microscopic physical operators appearing in H. These are 

ai = (F\ ai \F) = Y, V^(F l \Um-l) = J2 V^Tl F m r m+ n (4) 

(F\a ia +\F) = (FilailFi) (Fjla+lFj) = a,a*, and (F\(ni) s \F) = £ n n s \p n \ 2 where the 
exponent s is an integer. To calculate (\l/|c^|\l/) and 7i = (^\H\^) in (\E'|S'|\I / ) we 
standardly set = e zA \F), the phase A representing the effective action within the 
variational procedure. The first quantity becomes (\E f |c^| 1 I r ) = iA + (F\d t \F) with 



while H = (F\H\F) reads 

n = (E> 2 - ^i 2 ) - E«,> T ^ ^ ( 5 ) 

From the action A = j Ldt = j dt[ih{F\dt\F) — TC] (in which the second term of 
(F\d t \F), a total time-derivative, can be eliminated) one obtains the Lagrange equations 

d dL dL d dL dL 

didK~dE = 0, didK~dK = 0, 
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(df^/dt = /*J that can equivalently be written as (we set H = 1) 
By defining the Poisson brackets 

dA OB OB dA 



{A,B} = -iJ2 E 



(7) 



.0$ d/n r7/i! 

equations © can be formulated within the standard Hamiltonian formalism as dj l m jdt = 
{fmi'H} and df^/dt = {f^Ti.}. The resulting M.T dynamical equations are 

(equations for are obtained from the latter by complex conjugation), where 

and, in addition to definition HI one has a* = ^„ =0 y/m fmfL-v Concluding, we recall 
that within the |F)-based scheme the average total particle number TV = (F\N\F) = 
^,.(F|nj|F} = J2iJ2 n n \fi'\ 2 should be a constant of motion of Hamiltonian [5j The 
validity of this property is verified in Appendix C| In addition to A/", other M motion 



constants can be shown to be represented by F = |/™| 2 = (Fi\Fi). These allows one 
to implement the local-state normalization condition = 1. 

2.1. Connection with the Glauber-like trial state scenario 

A simple assumption allows one to recover state \Z) from \F) and relate the 
corresponding variational schemes. This is 

fin = e-WPzT/Vri (9) 

entailing Yln=o fn\ n )i = e _|Zl|2/2 J2n=o ( a tT\^)i/ n - arid therefore 

M oo M M 

\ p ) = nE*>* = n^ - ^! )* = n i*> = ^ > ^ 

i=l n=0 i=l i=l 

where the defining formula of Glauber CS \z) = e za+ ~ z * a \0) = e -\ z \ 2 / 2 e za+ |o) has been 
used for each space mode together with the Baker-Campbell-Hausdorf decomposition 
formula [1]. The same assumption enables one to find the new form of parameter H] 

oo oo I |2 m 

aj = E + 1 fUL+i = e ~ lZj]2 E z i ~~T = z i ' 

m=0 m=0 

showing that aij reduce to Glauber CS parameters Zj, and formulas (F\aia1-\F) = ZiZ*, 
(F\rii\F) = \zi\ 2 and (F\n 2 \F) = |zj| 4 + \zi\ 2 . In order to recover the M.T equations 
inherent in the |Z)-based picture we consider the time derivative of ctj. This is given 



Some remarks on the coherent- state variational 



6 



by iaij = Xlm=o V m + 1 [ifL+i df^/dt + if^ df % m+1 / 'dt], which reduces to (the detailed 



calculation is carried out in Appendix B) 

dcti 



m=0 



(12) 



Notice that, as illustrated in Appendix B , no explicit assumption on the form of has 
been requested so far (except for (Fj\Fj) = 1) in getting [12] At this point, however, 



the use of formula 191 in [T2l becomes necessary. We find ^^=o Umyrn + 1 



'm+l 



U z j Em=o m fmfm = Uzj\zj\ 2 which leads, in turn, to the well-known final equations 



IZi 



Uzj\ Zj \ 



(13) 



describing a set of discrete nonlinear Schrondinger equations [41], namely the M.T 
dynamical equations associated to the Bose-Hubbard model within the Glauber-like 
variational picture [6] . Equations [13] can be issued from the new Hamiltonian 

17 1 



n 



(14) 



{A,B} = -iJ2 f 



{zj, z}} = —idj 



(15) 



2 ^— ^^<^> 

obtained by rewriting formula [5] in terms of Zj, and defining the new Poisson brackets 
(PB) 

dAdB dBdA' 

dzg dze dz£ dz,£_ 

The crucial point that explains and justifies the whole reduction of the Hamiltonian 
picture based variables /£, /~, to the one involving the restricted set of variables Zi, Z£ 
thus consists in showing that PB[15]are consistent with PB of [3 In particular, a,, a} 
must be shown to exhibit, within algebra UJ the same algebraic structure of variables 
Zj, z\. By setting A = at and B = a\ in PB [7] one discovers that 

{a v a\} = £ E V^+TV^+l{fLf J m+1 ,fUn + i} = 



-is^ (m + 1) (i f j 2 - \f m+1 \ 2 ) = i/i 



(16) 



due to normalization = 1. Also, one easily proves that {aj, N^} = —idjiaj, where 

Ne = (F\ne\F) = J2 n n \fn\ 2 - Hence, it is an intrinsic feature of algebra [7] characterized 
by {fl,fm} — —ibjzbnm the property that aci, a} form a (classical) Weyl-Heisenberg 
sub-algebra of algebra [7J Noticeably, the latter represents the classical counterpart of 



the original boson algebra 



Sjeaj, 



l a j ) a £ 



Sj£ charaterizing Hamiltonian [3j 



Then identities a>j = Zj, and a\ = z\, obtained by assuming as a function of Zj (see 
formula [TT!) . quite naturally entail that zj, z\ obey the canonical brackets given in 



within the Glauber-like scheme. This completes the proof that the |Z)-based variational 
picture is consistently contained within that based on the more structured state \F). 

Concluding, we notice that if an effective Hamiltonian 7i depending on /£, can 
be rewritten in terms of collective variables atj, a} then sub-algebra [16] is sufficient for 
determining the evolution of the system and the |F)-based picture becomes redundant. 
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This is not the case of Hamiltonian where, owing to the presence of the nonlinear 
[/-dependent term, algebra [7] is necessary to derive the relevant motion equations. 

Comparing equations M and [13] fully evidences how the more pronounced quantum 
character of the |F)-based picture involves a dynamical scenario of great complexity. The 
marked semiclassical character of the |Z)-based picture instead appears when comparing 
quantum model [3] with Hamiltonian [TH The latter, in fact, is essentially obtained from 
[3] through the substitutions a { — > Zi, af — > Zi, namely by implementing the Bogoliubov 
approximation. At a formal level, the |Z)-based scheme thus provides an effective, 
dynamically-consistent formulation of the Bogoliubov semiclassical picture. 



3. Mean-field approach based on state |£) 

A quite significant form of state \Z) given by [10] is achieved with a simple calculation 

M 

is>=ni*> 

i=i 

5=0 b - \ i J 5=0 

where Il^ 1 |0) i = |0,0, ...0), & = Zi/y/N, and \S;£) corresponds to state |£) defined by[2] 
where the group- representation index S has been evidenced. State \S; £) is characterized 
by {£'■> S\N\S; £) = S and the orthogonality property ((; S'\S; t;) = 5gs>- The final 
version of \Z) derives from the observation that (Z\N\Z) = \zi\ 2 = M is the average 
total boson number in |Z)-based scheme and that £j = Zi/yfN is consistent with the 
normalization condition ^ |^| 2 = 1 of SU(M) coherent states. The latter follows from 
the scalar-product formula of two CS states |£) and \rj) given by (t]\£) = (EjT/*^) . 

The new information about \Z) is therefore that states \S;£) are its constitutive 
elements. In particular, state [2] features the property of incorporating only contributions 
of Fock states pertaining to the S'-particle sectors of the Hilbert space. This becomes 
quite evident rewriting \S',£) as 

\S;Z) = -7= iy^&<4\ |Q) = V -t =^ — er -Q M \rn) (17) 

where |0) = |0, 0...0) and superscript (S) recalls that S = Y^i m i anc ^ I 7 ™) ^ s sucn that 

\rh) = \mi, ...itim) = I I — ) 1 0, ...0) =>- N\m) = \ mArh) = S\m) . 

. \/mA ' / — ' 

The previous formulas allows one to evaluate the weight of state \L; in \Z) 
«; L\Z) = e'W ± ^ (C; L|S; = e Hf *J <•?,) ' . 

5=0 \ i / 



E, M : 



n 



|0)» 



E, 



kil 2 /2 e Ei 



'|0>i 
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Upon setting £ = £, the normalization condition ^\ |£j| 2 = 1 entails that (£;L|Z) = 
e ^ 2 Af L ^ 2 /\fL\ whose maximum value is reached for L = Af (Af is assumed to be 
integer). Considering \L;£) with L = Af ± p and p « L, one easily finds that the 
state-weight distribution around the maximum-weight state \Af; £) is not sharp. 

The variational procedure reviewed in section [2] can once more be applied to BH 
model [3] assuming |^) = e %A \£) as the trial state. The weak form of Schrodinger 
equation (^\(ihd t — H)\^f) = provides the action A = j dtC(£) where the effective 
Lagrangian £(£) = ih(£\dt\£) — (£\H\£) supplies the dynamical equations of variables 



The calculation of both (£|0t|£) and Ti.(i) = {£\H\£) is carried out in Appendix D 
together with the basic formulas required to achieve these results. We find, in particular, 
that the average local boson number is = Af\(,j\ 2 , reproducing consistently 

(£|iV|£) = N"Ylii\ij\ 2 = N. This suggests to define variables ipj = VN£,j (formally 
coinciding with Zj) for better comparing the present AAfF dynamics with that issued 
from trial state \Z). Explicitly, one finds (£|<9t|£) = NY^j£j€j = ^j^j^j an d 

3 <jt> 

while the dynamics is found to be governed by 

^ = ^^WV 3 -£r iA . (is) 

Equations [18] can be interpreted as the projection of equations [13] on a given S-particle 
Hilbert-space sector. In order to prove this property one has to consider the variational 
scheme based on a generic state \ip) = Y^s^s\S;C)- The latter reproduces state \Z) 
when the condition C$ = (S^^^e^^Af 3 ^ 2 is imposed. The |^)-based scheme would 
involve, in this case, the effective Lagrangian £ = ih(ifj\d t \ip) — (ip\H\ip) which, in the 
case \if)) = \Z), leads to equations [TBI Observing that ((,; R\S; (,) = 5rs, then 

wm) = "£j2 s c * c s& R \ x \ s > o = z2 s c s& s \ x ^ o > 

for both X = H and X = d t , states <9 4 |S;£) and H\S;£) pertaining to the S'-particle 
sector. Hence, C reduces to a summation C = J2s £s(0 over independent S'-particle 
Lagrangians £s(£) = (£; S\[ihd t — H]\S; £), the case S = Af giving equations [TBI 

Formally, no significant difference therefore distinguishes equations dH] (and the 
relevant generating Hamiltonian) from the picture corresponding to equations [13] if 
(Af — l)/Af — > 1 namely for boson numbers Af sufficiently large. A profound difference 
instead concerns the role of variables Zj and ipi in the relevant schemes. While 
Zi = (Z\di\Z) relates to the local superfluid parameter z i( its counterpart in the 
based scheme has no explicit physical meaning since (£|ai|£) = 0. State Oi\£) belongs 
in fact to the (Af — l)-particle Hilbert-space sector thus resulting othogonal to the Af- 
particle state With state \Z) this effect is avoided since \Z) is spread on the whole 
Hilbert space. The equivalence between the two schemes is restored in the case of two- 
particle operators z-iz* = {Z\aiOj\Z) being comparable with ifiif* = Af^iCj — (£,\ a i a f\0- 
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Of course case i — j describing local populations ($|nj|$), $ = Z, f is also included. 
Variables ipj = \ipj\e j thus acquire a physical meaning in terms of local populations 
= l^il 2 - The relevant phases 9j have no role unless one considers expectation 
values of operators a^af involving phase differences 6^ — 9j. 

3.1. Group-theoretic form of state |f) 

State [2] displays the particularly nice property to possess a fully symmetric structure 
mirroring the fact that all modes a m play an equal important role in defining This 
symmetry must be "broken" for proving that [2] has the standard group-theoretic CS 
form where |£) is generated by the group action on an extremal state (the choice of the 
latter entails the loss of the symmetric form). To show this our first step consists in 
proving that formula [2] can be rewritten as 

|f) = (Af\y 1/2 (A + ) Af \0, 0, 0) = (Af\)~ 1/2 E {at) M E + \0, 0, 0) (19) 

where E + = E' 1 and E is an element of SU(M) whose parametrization in terms of 
variables £j can be easily determined. The action of af on the zero-boson Fock state 
|0) := 1 0, 0, ...,0) is the standard one (af) p \0, 0, 0) = \frij\\..., 0, rij, ...) with = p 
while <2j |0) = 0. We point out that the choice of generating A + from a+ rather than 
from af is equally possible and simply entails choosing, in turn, one from M possible 
parametrizations for |f ) and the relevant form of E. This arbitrariness reflects the just 
mentioned symmetry of formula [2] For proving [19] we show that E a±E + = A + where 
E = e tS e tD is defined as 

At M 
S = > D = /2 ®i( a t a i + a t a l) > { D+ = D ) 



i=2 



with <pi e R and 9i G R. Upon setting J2l=2 ®1 := ® 2 -> standard calculations show that 



e lD afe %D = ^jli%' a l ( see Appendix E[ ) with y\ = cos 9 and = i9ksm9/9 if k ^ 1. 



A further action of e tS gives 



M M 

JS JD „+„-iD „-iS 

h 



aie ~e -=£ yj e^a+ = ]T ^ = A + , (20) 

3=1 3=1 

where fi = e^ 1 cos9, f& = i9ke l ^ k sm.9/9 and the action of factor e l ^ iTli in e tS = H i e i 'l >ini 
is described by e^^a^e - ^"* = e^aj. The identification A + = Ea±E + suggested by 
formula [20] is confirmed by the fact that the correct normalization of A + components 

& follows from Y,f=i 101 2 = cos2 6 + Ejl 2 (^/^) sin2 # = L Since a>i|0) = rn\0) = 
and thus S\0) = D\0) = 0, formula [T9l becomes 

10 = (M\r 1/2 E(at) M E + \0) = m- 1/2 E(ai)*\0) = E\Af,0,...,0), 

being E+\0) = e~ iD e- iS \0) = |0) and (at)"\0) = VM\Af,0, ...). Upon observing that 
E = e lS e lD = exp [e' s if e" !S ] e tS state |f) can be rewritten as 

|0 = e ^ lAr exp [e i5 iDe~ i5 ] |JV, 0, 0) = e^T{QW, 0, 0) , 
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where e iS ata%e lS = e % ^ 1 ^>ata\ entails that 



M 



e iS D e - lS = J2 (Q a +a £ + Ctafax) , T(() := e^= 



2 (CI <H a l +C,eaJa 1 ) 



with Q = 9ie l ^ e ^ >1 \ £ = 2, 3 M. Summarizing, we have found that 

|0 = -j={A + )«\U, 0, 0) = e**^T{CM, 0, 0) , (21) 

where T(Q is an element of group SU(M), which proves that stated up to an irrelevant 
phase factor, is generated by the group action of T((). The identification of T(() with an 



element of SU(M) is discussed in Appendix F By setting 4>i = 0, the relation between 
£i and (e is described by & = QsvuO/O whereas £i is fixed by |£i| 2 = 1 — ^2fL 2 \£,e\ 2 - An 
initial choice with A + = Ea^ n E + in formula [19] would have entailed T(() generated by 
Y^tLm (Q a m a e + (i a t a m) and the extremal state |0, ...J\f, ...) = (A/'!) _1/2 (a+) A/ '|0). 

As a final step, we prove that formula [21] is consistent with the group-theoretic 
definition of CS based on the notion of maximal isotropy subalgebra (MIS). Within 
the CS Theory [1] a class of CS is derived by identifying the (complex) MIS B of Q = 
su(M) and the related extremal vector IV'q)- The defining formula for |?/> ) states that 
a\ipo) = A |^o), A a G C, Va G Go where B D G := Go- The (complex) MIS naturally 
related to formula [21] is given by 

B = {h k ,ata k ,a+a k (k^£) : k,£ G [2, M]} , ([B,B}CB) 

whose generators are such that afa k \J\f, 0, 0) = a~^a k \J\f, 0, 0) = and generators 
h k form the Cartan (abelian) subalgebra. The vector satisfying the defining formula for 
|^o) is thus \J\f, 0, 0). According to the MIS scheme, coherent states are generated by 
the action on \tpo) of the elements of the quotient group G c /B where G c = expG and 
B = expB. The algebra that generates G c /B is in our case {a^ai : k G [2, M]}, which 
entails that a coherent state, up to a normalization factor A, has the form 

Ae Efc??fe ^ ai |Af,0,...,0) , e ^ ma * ai G G c /B . (22) 
State [21] precisely has this form. In order to check this property, one can observe that 
T(C) = exp \i Zli 2 (C k a+a k + C^a^] = e i9 ^ D + D+ ^, wit h D = Z k=2 C k a k /6 where 

CfcO r „+i ICfc' 2 



M M Aat/ . M 

M=EE^Mfl = E J ^r = 1 - 



fc=2 «=2 fc=2 

The exponent of T(() can thus be viewed as an element of su(2) in the two-boson 
(Schwinger) realization with generators J_ = a^D, J + = aiD + , J 3 = (D + D — a^a\)/2 
and commutators [J3, J±\ = ± J± and [J + , J_] = 2J 3 . This information allows us to 
apply the standard decomposition formula e vJ +~ v * J - = e uJ + e in(i+|ix| 2 ) e -u*J- f or ^ ne 
SU(2) elements [1] where i>, w G C, v — \v\e 1 ^, u = \u\e 1 ^ and \u\ = tg\v\. Setting v = i9, 
which entails u = itg8, one has T(() = e %e{ - Dl B+D+a ^) = e iO(J-+J+) thus obtaining 

T(C)|A/ "' °' 0) = (l + \u\^ °' 0) = (1 + \ u \^ 0, 0) 
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where rj k = u( k /6 = ie i6k tg6, and J 3 \Af, 0, 0) = (7V/2)|7V, 0, 0) has been used 
together with J-\Af, 0, 0) = a^D\Af, 0, 0) = 0. Therefore state I2T1 indeed can be 
cast into the CS form 1221 determined within the theory of CS. 



4. The duality property of states \Z) and |£) 

Both state \Z) and whose definition involves boson operators aj and Oj of the 
ambient lattice, can be shown exhibiting a dual character which becomes evident when 
space-like operators are expressed as momentum-like operators through Fourier formulas 

^ e -iqj ^ e iqj 

h « = H^m a ^ ^ = £-7^' q:=2nq/M,qE[l,M} (23) 
j=i V j =1 V 

where [a,-, a/] = 5je implies that [b q , bp] = S qp . Notice that we assume periodic boundary 
conditions (namely the lattice is a closed ring) so that displacements q q + sM and 
j —* j + rM (r, s are integer) leave operators aj and b q unchanged. This condition, 
standardly assumed to simplify theoretical models, becomes necessary in real lattices 
with a ring geometry [42] . Concerning state \Z) = JJ. \zj) simple calculations yield 

\Z) = e-lE.fePgE^^+io) = e -|E fc KI 2 e E fc ^^| ) = JJ \ Vk ) = \v) (24) 

k 

where a^O) = = 6fc|0) has been used (recall that |0) = |0,0,...,0)) and, thanks 
to definitions [231 one has v k = J2jLi e"" 1 ^ Zj / \/M , and Zj = Ylf=i e ^ 3 v q /\fM. Trial 
states \Z) are thus equivalent to states \V) formed by momentum-like Glauber CS 
\v k ) = e Vkbk ~ v k bk \0) k . Similarly, states |£) transform into momentum-like SU(M) CS 

l ^' = ~7]k lo) = ~7W lo) = lAf,a) ' ^ = ^7W ak (25) 

where the latter definition ensures B + = J2h=i ak H = Sj=i£j a j~ = ^ + ■ Also, the 
counterpart of formula [T7] in the momentum picture is easily derived 

(AO 

\M;a)=J2 C ^<...a^\p), 

p 

where C^N) := [Af\/UkPk ] -} 1/2 while \p) = [UkPk ] -V 1/2 Uk( b k) Pk \ ) are momentum 
Fock states. While space-like states \Z) and |£) provide information on the local 
boson population by means of {Z\a^ai\Z) = \zi\ 2 and (£\afcii\£) = A/"|^| 2 , respectively, 
momentum-like states | V) and \a) provide information on the fc-mode boson population 
by means of (V\b^bk\V) = \vk\ 2 and (a\b^bk\a) = Af\ak\ 2 , respectively. The total boson 
number Af is unchanged being (Z\N\Z) = (V\N\V) and (£\N\£) = (a\N\a). 

As an application of the duality property of we show that states \Sk), describing 
Schrodinger cats, can be defined having specific momentum properties. To this end an 
important preliminary condition consists in showing that A/"-boson states \£(£)) with 
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£ G [1, M] can be exploited quite easily to form sets of M orthogonal states. Recalling 
that the scalar product of two SU(M) CS is given by (r}\£) = (YljVj^j) ■> one nas 

(mm) = (Eot&w)^ = s m Eseo&o = 

which displays how the desired orthogonality directly ensues from the orthogonality of 
complex vectors £(£) = (6(^),6(^)> •••) witn ^ e [1, M]. For fully localized states 
characterized by = the orthogonality condition (£(h)\£(£)) = She is manifest. 
In the general case, however, the condition = $h£ with \£e(£)\ >> \£j(£)\ 

can be achieved by exploiting the arbitrariness of the phases of £j(£). States | 
describing strong boson localization have been employed to realize Schrodiger-cat states 
\S k ) in a ring of attractive bosons [25J. These were proven to well approximate the low- 
energy states including the ground state in the regime of strong interaction. Following 
the recipe given in |25j we define \S k ) as a superposition of equal- weight localized states 

M iki 

e=i v M 

As a consequence of the orthogonality of states \^{£)) states \S k ) appear themselves to be 
orthogonal namely (S q \S k ) = 5 qk . We observe that, if (£(•£) Kl£C0) = A%(f)| 2 ~ M8a 
evidences the information about boson localization at the £th site, the expectation value 

M M i(%i-qh) kf 

(s k \ ni \s k ) = e E MQMtw = -h ' v? ' 

h=l 1=1 

obtained through the properties (^a^a^) = NVm&ivlO 1 ' 1 ^ an d (£WI£C0) = <W> 
confirms the the expected feature of full derealization typical of Schroedinger states. We 
note how the possibility of constructing set of orthogonal trial states is quite important 
for applications to boson lattice systems such as model [31 While trial states can be 
used for approximating in a reliable way sets of energy eigenstates, the possibility to 
make them mutually orthogonal certainly enriches the approximation with an important 
feature. 

In order to show that states \S k ) have specific momentum properties we rewrite 
in its dual form \a(£)). Thanks to formula [25] ]£(£)) = \a{£)} with a k (£) = 
Ef=x e-W&eyVM, one obtains = (p\a(£)} = C^M) < ...a™ giving 

(p\S q ) = -±= E e m (Pim) = 4^ E ^WN) a?(£) ...a&(£) . 



In case of strong localization condition \(,e(£) \ — 1 >> |£jCOI leads to the approximation 
a k (£) ~ e~ M j\[M and, in particular, to 

(fAc\ 1 S^ P Wr; (\r\ e _ V - /V! r -^y,An( Pk \) J^\q-Mp)] 
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where C^N) := [AA!/n fe p fc !] 1/2 and \{p) = ^ k k Pk = 0,1,2, ... mod(M). The latter 
represents the eigenvalue of the total quasi-momentum operator P = a k b\b k with 
a = 2n/M such that P\p) = o~\(p)\p). It is worth recalling that, in the discrete 
geometry of ring lattices, the quasi-momentum properties are described through the 
displacement operator D = exp[— iaP), whose action is displayed by DaiD + = a^+i and 
Db k D + = b k e lkcj . Based on equation D\p) = e~ taX ^\p) Fock states can be organized in 
M equivalence classes labeled by X(p). Index q in \S q ) therefore characterizes the quasi- 
momentum associated with \S q ) since the term J^exp [i£[q — A(p)]] in (p\S g ) vanishes 
whenever \p) has a momentum X(p) ^ q mod(M). States \S q ), in the presence of 
strong localization, supply a set of M orthogonal states whose label q bears information 
on the class with quasi-momentum A maximally involved in the realization of \S q ), 
within the TV-particle Hilbert space. As for Glauber-like states, we note that (X\Z) = 
Ylj( x j\ z j) = Ylj ex P[^j z j ~ (\ z j\ 2 + l x il 2 )/2] so that \Z) and \X) cannot be orthogonal. 
At most, M quasi-orthogonal states can be obtained by considering sets {xj(£)} such 
that Ix^f)! 2 ~ M8jt for which \ (X(h)\X(£)) \ ~ e~ N ' . Representation of Schrodinger-cat 
states in terms of quasi-orthogonal states \Z) can be developed under these conditions. 

5. Conclusions 

In this paper we have compared variational schemes based on trial states \F), \Z) and 
widely used in applications to many-mode boson systems. To illustrate their distinctive 
features we have applied such schemes to the BH model which has became, in the 
recent years, the paradigm of real interacting-boson lattice systems. Such a comparison 
has been aimed at evidencing the specific characters of each scheme to favour their 
applications in the study of the properties of many-mode boson systems. 

In section [2, we have applied the |F)-based scheme to the BH model showing, 
within the corresponding dynamical scenario, that collective variables Oj, a* form a 
(classical) Weyl-Heisenberg sub-algebra in the Poisson algebra of variables fl, f^. This 
crucial property allows reduction of the |F)-based picture, exhibiting a more pronounced 
quantum character, to the |Z)-based picture based on Glauber's CS. The semiclassical 
character of the latter appears to be an effective, dynamically-consistent procedure 
incorporating the Bogoliubov approximation. 

In section [3j we have shown that Glauber-like trial state \Z) is a superposition 
of SU(M) CS |jV, £) that involves all the A/"-particle sectors of the Hilbert space. We 
have exploited this property to explain why the dynamical equations relevant to the 
BH model obtained in the \J\f, £)-based scheme coincide with the equations derived in 
|Z)-based scheme. The meaning of microscopic CS parameters of such schemes has been 
illustrated and related to the fact that states \Af, £) are boson-number preserving. Also, 
in section [31 we have discussed explicitly the procedure that enables one to recast state 
[2] into the standard form |£) = g\Q) of CS theory involving the extremal vector 

Section H] has been devoted to discuss the duality property of space-like states 
\Z) and |£) which allows to rewrite them as momentum-like states involving modes 
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bt. This property has been used for constructing Schrodinger-cat states with specific 
momentum features for bosons in ring lattices. In general, the use of states |£) and 
the possibility to construct set of orthogonal states outlined in section @] should allow 
a better characterization of low-energy regimes in systems of bosons confined in ring 
lattices whose standard description is given in terms of Hamiltonian [3j Particularly, 
the duality property of states \Z) and |£) finds a natural application in the study of 
supercurrents and vortex states occurring in such systems [22], [23J. 

As observed in section (TJ while the |Z)-based scheme has been extensively used 
in applications, the interest for the |£) and |F)-based schemes is more recent. The 
more pronounced quantum character of the |F)-based scheme is expected to supply, 
in the applications to BH-like models, a better description of the critical behaviors 
[29]- [32] inherent in quantum phase transitions. For the same reason it should supply 
an effective tool to study the complex dynamics [27], [28] of bosons in lattice systems. 
Bosons distributed in small arrays (and thus involving small number M of space modes) 
are specially interesting since they can switch from fully quantum to intermediate 
semiclassical behaviors by adjusting model parameters [23j-[25j. The correspondingly 
small number of components \Fi) in \F) makes it feasible for performing numerical 
simulations of equations [SJ These aspects will be investigated in a separate paper [43J. 

Appendix A. Application of formula [2] to the case M=2 

It is quite easy to show that formula [2] with M = 2 reproduces the standard group- 
theoretic definition of SU(2) coherent state. In this case A + = + £,2 a t- Then 

W > = (A/1)- 1/2 (A+nO, 0) = g (4)>+)^|0, 0) 

= E C.W -8,s) = e"« J2 C 'W* | J; - J + s) = e"*\ J; z) 

s=0 s=0 V "i" I N 

where C s (J\f) = VAf\/ \J s\{M — s)\, the definition z = nas been used, and <\>i is 

the phase of Moreover, | J; — J + s) = \Af — s, s), where J = Af/2, can be seen as the 
mth vector (with m = — J + s) in the standard basis {| J;m) : J^\J\m) = m\ J; m)} of 
algebra su(2) within the Schwinger boson picture of spin operators J3 = (a^~a 2 — a^"ai)/2, 
J + = G + fll = (J_)+. State g\Q) = e^ J +^ J -\J; -J) = (1 + \z\ 2 )~ J e zJ+ \J; -J), obtained 
through the standard decomposition [2] g = e^ J+_ ^* J_ = e zJ+ e J 3M 1 +\ z \ 2 ) e ~ z ' : J- ) coincides 
with state | J; z) just defined, where z and ( have the same phase and \(\ = tg\z\. 



Appendix B. Derivation of dynamical equations for Zj 

Based on dynamical equations [8] governing the evolution of variables one has 



.dot 



dt t-^m. 



,d_l 
dt 



m rt 
J 1 



m+1 



+ i 



df, 



m+l fi 



dt 



J n 



Some remarks on the coherent- state variational 



15 



V UmVm + lfUL+1 + ®i \ J2 (i m + l)/i+l/m+l - ™ /m/m) ~ E /m/S 



+ 



/"m+2 



where the index- m range is [0, oo]. Thus dj is formed by three terms. Substitution m — > 
m+ 1 in the first summation of the third term (one should notice that = 0) shows that 
the latter vanishes, while, in the second term, S m [(m + l)/m+i/m+i ~ m fmfm\ = 
is easily proven. Further simplification is achieved if the the on-site normalization 
condition (F^Fi) = J2 m fmfm = 1 is imposed. Notice that choice M automatically 
ensures such a condition since = \zi) and Glauber CS are such that {Zi\Zj) = 1. 
Under this circumstance the second term reduces to 



Appendix C. Conservation of TV and other constants of motion 

After recalling that M = (*|JV|*) = (^\J2j n j\^) = EjEn™!/™! 2 where = 
Ylm n \fn\ 2 an d n e [0' °°]' i e [1) ^ us consider the Poisson bracket of M and H 

w«} = E-E *{i#i a ,«} = EE = 

• /;:' - ")/;!' + v^+T^+i ^ + vW,! ; -i-;) 
= -*E/;E„ (^ATFT/^/i) -*E^E n (-v^n^ + i^) 

= -* E; I> i n ' - E/^* = ~* E^E tej ( T « a < - = • 

Then dM/dt = {Af,H} = 0. A similar calculation allows one to evidence that other 
M constants of motion are involved in the dynamics of /•[. These are Ij = En I /n 1 2 f° r 
which {Ij, 7i] = 0. Quantities Ij are in involution with A/", and ij namely {ij, J\f} = 
and {/,-, Oij} = Vj. One easily check as well that {/_,-, 1} = for each i and j. 



Appendix D. Formulas relevant to the SU(M)-CS picture 

In order to obtain £(£) = ih(£\dt\£) — ((,\H\t;) one needs to calculate (£|<9t|£) and 
(eKni-lKIO mH(0 = {£\H\£). Concerning ^|(^-l)n^), (recall that |f> = \S;0) 
one should observe that [oj, (^4 + ) s ] = s£j (A + ) s_1 , and that 

O.I0 = a lPs (A + ) s \0) = Ps [(A + ) s a t + ^(A+) 5 - 1 ]|0) = Vs&O 
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where ps = l/V~Sl, and = ps_i(A + ) 5_1 |0) is a (S*— l)-boson coherent state. Iterating 
this calculation gives af|£) = gy/S(S - 1) with = p 5 _ 2 (^ + ) 5 ~ 2 |0). Therefore 
(CK|0 = ... = iV|&| 2 whereas 

(CKn, - l)n»|0 = (CIK + )X 2 |0 = S(S - 1) |^| 4 (£"10 = S(S - 1) |&| 4 . 

Similarly, one finds (^a+a^) = (0\p s A s aia tPs (A+) s \0) = S^M = For 

two generic states \r)), |£) the latter becomes (^a^a^) = Sr/^i (£V ?7*£i) s '~ 1 where the 
inner product (r)\£) = (Xli^*^) 5 °f two 'S'-boson states |£) has been used. In the 
effective Lagrangian C term (£|9 t |£) can be recast as 

m\o = = pLi(oi^ 1 E,^<(^ + )^ 1 i°) = «i£^;io 

= <^E^ + s] io = E.,^1 + E,E m c^ w^io 
= £>s + E,E m c4 (s - ous = - - oE.&s = 5 E^l • 

Concluding, the four /two-boson expectation values just obtained provide the effective 
Hamiltonian (£|if|0 = f S(S - 1) £\ |^| 4 - S occuring in £(£)■ 

Appendix E. Derivation of operator A + 

After setting ®1 := ^ the action of e lD on a| is given by the standard formula 

oo . s 
s=0 S ' 

where [D,a+] s = [D, [D, a+] s _i], [D, af]i = [_D,a+], and [D,a+] = 1. Observing that 
[D, a^] 2r = 2r ai" and [D, af] 2r +i = 9 2r Q with Q = J2h=2 $fc a fc > one obtains 

oo (_\r °° •/ \ r • /] M 

e-a+e-- = £ j-U V + E T^TlV ^ = ^ cos 5 + = £ 

where yi = cos 6* and y& = i6 k sm6/6 if A; ^ 1. 
Appendix F. Two-boson operators of algebra su(M) 

The fact that T(£) G SU(M) is easily demonstrated by recalling that, within a Schwinger- 
like picture, algebra su(M) can be realized in terms of two-boson operators a^a^, a^aj 
with 1 < j < M — 1 and j + 1 < k < M that play the role of lowering and raising 
operators, respectively. This set is completed by the generators of the Cartan-subalgebra 
{hk, k — 2, ...M : [h k , he] = 0} where each of the M — 1 operators /i^ can be written as 
an appropriate linear combinations of number operators Hi = afai, % = 1,2... , M. We 
notice that, consistent with the presence of the group-invariant operator N = n ii 
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only M — 1 operators h k can be realized with M operators ?v A generic element of Q 
= su(M) = {a+a m (m ^ j) : m,j G [1, M]; : G [2, M]} is thus given by 

M-l M M 

Yl Yl ( Z kJ a t a k + Z lj a 3 a t) + Y akhk 
j=l k=j+l k=2 

where z k j = x k j + iy k j and x k j, y k j, a k G R. Since elements g G G of a Lie group 67 
are generated by the Lie algebra element a G Q = Lie(G) through the exponential map 
g = exp(ia) then the latter fromula shows that T(t) G SU(M). 
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